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Studies of angle-dependent magnetoresistance oscillations (AMRO) in the interlayer conductivity 
of layered metals have generally considered semi-classical electron transport. We consider a quantum 
correction to the semi-classical conductivity that arises from what can be described as an interlayer 
Cooperon. This depends on both the disorder potential within a layer and the correlations of 
the disorder potential between layers. We compare our results with existing experimental data on 
organic charge transfer salts that are not explained within the standard semi-classical transport 
picture. In particular, our results may be applicable to effects that have been seen when the applied 
magnetic field is almost parallel to the conducting layers. We predict the presence of a peak in the 
resistivity as the field direction approaches the plane of the layers. The peak can occur even when 
there is weakly incoherent transport between layers. 

PACS numbers: 71.18.-|-y, 72.10.-d, 74.72.-h, 74.70. Kn 
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I. INTRODUCTION 

Angle dependent magnetoresistance oscillations 
(AMRO) are a valuable tool for the investigation of 
properties of layered metals such as organic and cuprate 
superconductors. The dependence of the interlayer mag- 
netoresistance on the direction of the magnetic field has 
been used to map out a three dimensional Fermi surface 
sometimes in exquisite detail in a diverse range of layered 
metals. These include organic charge transfer saltsyi 
strontium ruthenateSf^s^ semiconductor superlattices?^ 
monophosphate tungsten bronze,^ and an overdoped 
cuprate-^— More recent developments have extended 
this to allow one to not only determine the anisotropics 
in the Fermi surface, but also of the anisotropy of 
scattering on the Fermi surface ^^i^iiSiiSiiiiiHiH 

In addition to the shape of the Fermi surface, the 
coherence of interlayer transport in layered metals has 
been a complicated and controversial issue. ^"^ In order for 
AMRO to be observable, only coherence between neigh- 
bouring layers is required, i.e. weakly incoherent trans- 
port, and for almost all angles, this leads to AMRO that 
are equivalent to fully coherent transport (i.e. there is 
a three dimensional Fermi surface) perpendicular to the 
planesjiiii^ In general the two situations can be distin- 
guished in transport measurements only if there is a co- 
herence peak in the resistivity for fields close to parallel 
to the layers. In this work, we show that quantum inter- 
ference effects can also lead to a peak in the resistivity at 
fields close to parallel to the layers, for weakly incoherent 
transport. However, this peak can be distinguished from 
the coherence peak by its magnetic field and temperature 
dependence, which is determined by the lengthscales and 
timescales over which quantum coherence is destroyed. 

The peak in the resistivity for parallel fields that we 
find arises from a contribution to interlayer conductiv- 
ity in a magnetic field B = i3(sin 9 cos (j>, sin 6 sin (p, cos 9) 



in layered metals which is a quantum correction to ex- 
isting semi-classical transport formulae.— i^^ There has 
been previous work in this direction j^^i-'^'^'^^i^^'^^i^^'^^ but 
this was either focused on weak localization in anisotropic 
metals, or did not consider the magnetic field configura- 
tions relevant for AMRO. Our work is motivated by ex- 
periments in which deviations from conventional AMRO 
are seeUi ^^i^'^i^^i^^i^^ The particular deviations of interest 
are weak-localization like peaks at angles near 9 = 90° j^ 
We calculate the quantum contribution to the interlayer 
resistivity and show that it behaves in a similar way in 
a perpendicular magnetic field to the Cooperon in a two 
dimensional electron system. We also show that a peak 
for parallel magnetic fields emerges naturally for reason- 
able choices of relevant materials properties, leading us 
to suggest that this physics is relevant for a number of 
recent experiments. 

Our main results are an expression for the inter- 
layer Cooperon correction to interlayer conductivity, and 
asymptotic expressions for magnetic fields that are close 
to parallel to the conducting layers. We demonstrate 
the presence of a small peak in the resistivity for such 
fields, and discuss the characteristic behaviour with tem- 
perature and perpendicular magnetic field that allow this 
feature to be distinguished from coherence peaks at sim- 
ilar field orientations. 

We note that the calculations we consider here may 
also be of interest to those considering weak localization 
effects in bilayer or multilayer graphenC j^^'^^ semicon- 



ductor double quantum wells )« semiconductor quantum 
wells with intersubband scatteringf^^ and superlatticesi^ 
The paper is structured as follows: in Sec.|TT]we briefly 
review the physics associated with magnetoresistance due 
to weak localisation and discuss the experiments that mo- 
tivate our calculations. In Sec. lIIII we give details of our 
calculations and results, and in Sec. IIVI we discuss the 
consequences of our results for interpretation of AMRO 



experiments, in particular on the issue of coherent versus 
incoherent interlayer transport. 



II. MAGNETORESISTANCE DUE TO WEAK 
LOCALISATION 

In order to connect our calculations to the experimen- 
tal data showing unusual peaks in AMRO, we first re- 
view the key physical quantities which determine the 
magnitude of weak localization effects in two-dimensional 
metals.— 1^ The total conductivity of the system in a per- 
pendicular magnetic field B± can be written as 



a{B^,T) = a,, {B^ , T) + 5a{B^ , T) , 



(1) 



where the first term is the semi-classical contribution and 
the second term is associated with quantum interference 
effects and weak localization. The relative magnitude of 
the two terms is of order 



(2) 
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where kp is the Fermi wave vector and I — vpT i& the 
elastic mean free path, with vp the Fermi velocity and r 
the elastic scattering time. Hence, quantum effects will 
be enhanced by increased disorder, which reduces kpL 
At low temperatures in zero magnetic field 
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(3) 



where r^ is the inelastic scattering time, which includes 
all inelastic processes, e.g. electron-electron and electron- 
phonon scattering. Since generally, l/r^ ~ T^, with p 
some positive number, the resistivity diverges logarith- 
mically with decreasing temperature. 

The magnetic field scale (the phase breaking field) at 
which quantum interference effects are destroyed is B^ — 
h/{AeDT^), where D = ^VpT is the diffusion constant. 
The area associated with one quantum of magnetic flux 
is related to the Thouless length, Lth = \/Dt^. This 
magnetic field scale can be written in the form 

in order to emphasize the fact that this field scale is tem- 
perature dependent. Destruction of the quantum inter- 
ference effects by the magnetic field leads to a negative 
magnetoresistance which is quadratic in the magnetic 



field for small fields ,^i^i22i33 ^j^g^ d ^ n_ ^ ^ 
and T^ ':^ T (Bq ^ B^{T)) the change in the quantum 
contribution to the conductivity due to the magnetic field 
can be written in terms of ipiz), the digamma function. 



Act = a{B±_) - a{B±_ = 0) 
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(5) 



where the last equality only holds when B(f,{T) ^ B^. 
Equation ([5]) can be used to fit experimental data for 
the magnetoconductivity with one free parameter, Tip{T), 
provided the diffusion constant is known. For thin metal- 
lic films and semiconductor heterostructures this has 
proven to be a powerful method for determining the ab- 
solute value and temperature dependence of the inelastic 
scattering ratei^i^ 



A. Summary of Results 

Our calculations establish the following: 

1. When the temperature is low enough that the inelas- 
tic scattering rate is much less than the elastic scattering 
rate, i.e. r^ ^ r, with no magnetic field perpendicular 
to the layers present, there will be a small increase in the 
interlayer resistivity which is approximately logarithmic 
in temperature. 

2. At low temperatures there will be a small negative 
magnetoresistance on the scale of perpendicular magnetic 
fields of order the field Ba ^ B^, (defined in Eq. pTjl ). 

3. Even when the field parallel to the layers is large there 
will still be a contribution to the resistance due to weak 
localization similar to that which occurs in zero field. As 
the field is tilted at an angle 69, slightly away from the 
layers, there is a small perpendicular field, B± ~ BdO. 
Consequently, a graph of the interlayer resistance versus 
9 will then have a narrow peak at 9 — 7r/2, which has 
a width proportional to Ba/B when r^ > r. The tem- 
perature dependence of Ba implies the width of the peak 
is also temperature dependent. The height of the peak 
relative to the background magnetoresistance will only 
depend weakly on the strength of the field parallel to the 
layers for fields larger than Bp (defined in Eq. ([5])). 

4. The resistance peak is quite distinct from the "co- 
herence" peak which can occur ai 9 = ^ if there is 
coherent interlayer transport (i.e., a three-dimensional 
Fermi surface) j^i^ This coherence peak only exists at 



very high fields (i3|| 3> 'm*kpc^/ep/t±/eT{T)), has a 
width 69 = 2kpc{t±/ep) which is independent of temper- 
ature and field, and its height relative to the background 
scales with B and Tr^j^ 

5. There are two regimes of parallel field, separated by 
an intermediate field scale. 



^" = ^1- 



(6) 



For _B|| <C Bp, the corrections to the conductivity are 
essentially independent of i?|| , whereas the semiclassical 
background does depend on i?||. For i?|| > Bp, the cor- 
rection to the conductivity has the same B\\ dependence 
as the semiclassical contribution, so that the relative re- 
sistivity correction (i.e. the shape of the peak) becomes 
a function of B± only near 9 = ^. 



B. Brief review of experimental results 

Below we discuss several layered materials in which 
peaks in the angular dependence of the interlayer magne- 
toresistance have been seen for parallel fields, and discuss 
whether the evidence points towards weak localization 
peaks or coherence peaks. 



1. l3"-(BEDT-TTF)AAM(C20i)iY 

One class of layered metals where weak localization- 
like peaks have been observed^ is the family /?"- 
(BEDT-TTF)4AM(C204)3Y of organic charge transfer 
salts, where Y is a solvent molecule, M=Ga,Cr,Fe, and 
A^HaO"'". These materials have significant structural 
disorder, especially in the anion layers. ^^ This can been 
seen from the fact that the ratio of the resistivity at room 
temperature to that at low temperatures is of order one. 
In contrast, in many metallic organic charge transfer salts 
this ratio is as large as one thousand)^ 

The Shubnikov-de Haas (SdH) oscillations yield an os- 
cillation with period 330 T, corresponding to a circular 
Fermi surface with wave vector kp — 1 nm~ ^ . The cor- 
responding Fermi surface area is an order of magnitude 
smaller than in many ET materials. The small Fermi sur- 
face area arises because the fi" crystal structure leads to 
four bands, three of which are almost completely filledri 
so the Fermi energy is near the band edge. The effective 
mass associated with these SdH oscillations is close to 
one free electron mass. The corresponding Fermi energy 
is about 35 meV. The Dingle temperatures, which are a 
rough measure of the disorder in the sample, estimated 
from the SdH oscillations are in the range 1-4 K, up to 
an order of magnitude larger than the cleanest organic 
charge transfer saltsji 

These materials have several features unlike other or- 
ganic charge transfer salts that make them more likely 
to exhibit weak localization effects: (a) smaller Fermi 
wavevector, (b) stronger disorder, (both (a) and (b) re- 
duce /cf^), and (c) the almost 3/4 filling will enhance the 
effect of the strong electronic correlations which narrow 
the electronic bands, reducing the Fermi energy e^?,?? 

In both the M=Ga and Cr (Y=C2H2Cl2) materials 
in zero field the interlayer resistivity increases approx- 
imately logarithmically with temperature below about 
20 Ki^ The total change is about one per cent in the 
temperature range 1-20 K. At a temperature of 1.5 K, 
perpendicular fields of about 0.2 T and 2 T, destroy this 
feature in the M=Ga and M=Cr materials respectively. 
Furthermore, as the field is tilted towards the plane of the 
layers there is a small peak sX 9 = 7r/2 in the interlayer 
resistance versus angle for both the M=Cr and M=Ga 
material. The width of this peak, for the M=Cr ma- 
terial, 59 (in radians), has the field dependence B^/B, 
with B^ = 0.25 T, between about 3 and 14 T^ Note 
that this is the same scale on which the weak localisation 
peak occurs for 9 = and small fields. The magnitude 



of the peak is about 0.3 per cent of the background mag- 
netoresistance and depends weakly on the parallel field. 
As the temperature increases from 0.7 K to 3 K the peak 
becomes broader and smaller, and is not visible at 6 K, 
whereas the background magnetoresistance changes little 
in the same temperature range. 



2. a-[(ET)i.^(BETS)^]2KHg(SCN)4 

In a-[(ET)i_:,(BETS):,]2KHg(SCN)4 with x ~ 0.03, 
the resistance versus temperature curve shows an upturn 
below about 5 K.''*' About ten per cent of this growth 
can be suppressed with a magnetic field. The peak that 
occurs at 9 = tt/2 is about 50 times broader than the 
coherence peak that occurs in clean a-(ET)2KHg(SCN)4 
samples at fields of order a few tesla<^ For low fields, up 
to 1.5 T, the shape of the peak depends only on B±, as 
we would expect for weak localization behaviour, and at 
higher fields, the peak becomes broader than at low fields. 
The high field behaviour that is not captured in our the- 
ory may be related to approximations we use that should 
break down at very large fields (e.g., 15 T in this case). 
In contrast, the coherence peak for a-(ET)2KHg(SCN)4 
has an absolute height that increases monotonically with 
increasing magnetic fieldi^i^ 



3. K-(BEDT-TTF)2Cu(NCS)2 

The magnitude and temperature dependence of the in- 
elastic scattering rate was recently determined in this 
material using AMRQji^ The data was fit to the form 



nT) Til 



(7) 



with Til — 3 ps and A ~ 0.006 ps^^K^^. This temper- 
ature dependence is also consistent with that of the dc 
resistivity in many organic charge transfer salts. "^^ The 
first term is associated with elastic scattering due to dis- 
order. The quadratic temperature dependence can be as- 
sociated with inelastic scattering due to electron-electron 
interactions in a Fermi liquid. The ratio of the inelastic 
scattering time t^(T) to the intralayer elastic scattering 
time is t^/tu = 50(K/T)^. Thus, as the temperature 
increases from 1 K to 20 K the ratio decreases from 50 
to 0.1 and Tcj,/tii will decrease from about 1 to 0.1. 

We also suggest that weak localisation may be the ori- 
gin of the feature near 9 = ^ that Singleton et al}— 
assigned to the coherence peak. The peak width appears 
to broaden with increasing temperature (as would be ex- 
pected for decreasing r^, and hence increasing B^(T)). 
In earlier data on the same material, ^"^ the peak height 
drops with increasing temperature and the peak width 
increases with increasing temperature. For an inelastic 
scattering rate that increases with temperature this is to 
be expected for a weak localisation peak, since B^{T) 



increases with increasing temperature and the width of 
the peak goes as B^, whilst its height should increase 
logarithmically with temperature. 



4- Intercalated graphite 

A peak near parallel field has also been seen in interca- 
lated graphite materials which display unusual inverted 
AMR04^ However, the width of the peak appears to be 
independent of magnetic field which would tend to argue 
against weak localization effects. 



the AMRO from weakly incoherent interlayer transport 
could be derived using a "Diffuson" -like equationiii We 
expect that there should also be a contribution to the 
inter-layer conductivity from a "Cooperon" like process, 
that will have a different magnetic field dependence to the 
"Diffuson" term due to their differing time-reversal prop- 
erties. It is this question that we address in this paper, 
and we find that indeed, an extra quantum correction to 
the interlayer conductivity is generated. Moreover, it has 
distinct experimental signatures. 



B. Impurity correlations in a layered metal 



III. MODEL AND CALCULATIONS 

A. Interlayer charge transport 

In previous work^ we considered the AMRO that 
arise when there is an anisotropic Fermi surface and 
anisotropic scattering, in two different limits. These two 
limits were the limit of coherent interlayer transport and 
the limit of weakly incoherent inter-layer transport in 
which there is hopping between adjacent layers. We shall 
not consider anisotropy in the Fermi surface, scattering 
or the interlayer hopping here, as these are additional 
complications beyond the physics that is our main inter- 
est. For isotropic weakly incoherent interlayer transport 
the Hamiltonian describing interlayer charge transport is 



n = t^ 



E [4^ 



h.c 



(8) 



where t± is the interlayer hopping integral. We assume 
that we are in weakly incoherent regime such that 



1 

ti_ < -. 
r 



(9) 



There is no observable difference in the AMRO between 
the coherent and weakly incoherent models of interlayer 
transport except for polar angles of the magnetic field 
very close to 90°, corresponding to a magnetic field par- 
allel to the planes.— J^- The interlayer conductivity that 
is deduced from such a hopping Hamiltonian isr^ 






+Gii(ri,r2)Gf(r2,ri)), (10) 



where c is the interlayer spacing, L is the sample size, 
and Gj is the retarded (advanced) Green's function 
in layer i. The angle brackets (. . .) indicate an average 
over disorder. Generically in calculating the disorder- 
averaged conductivity, there are two classes of diagrams 
to consider, the ladder diagrams that lead to the "Diffu- 
son", and the maximally crossed diagrams that lead to 
the "Cooperon" ,^2i^ In our previous work, we found that 



We now discuss the properties of the probability dis- 
tribution for impurities in a layered metal, and their con- 
sequences for impurity correlations. We assume that the 
probability distribution p(r) for the locations of impuri- 
ties can be factorized into separate distributions describ- 
ing the distributions parallel (pn (x, y)) and perpendicular 
{pz{z)) to the layers2i 

p(r) =p||(a;,j/)pz(z). 

We assume that an individual impurity (or lattice defect) 
located at R has a potential Vimp (r — R) associated with 
it, and that 



(^(r)) 



dVp(r')Mmp(r-r') = 0. 



Correlations in the random potential take the form 



{V{r)V{r')) = n, 



d3rXr")T^„,p(r-r")M,„p(r'-r") 



n,; / £v"p\\{x",y") j dz"p,{z") 
xl^,„p(r-r")yimp(r'~r"). 



(11) 



We shall assume that there is a uniform distribution of 
impurities in the planes parallel to the conducting layers, 
i.e. p||(x,2/) = -p-, where L is the sample size. All struc- 
ture in the impurity distribution is thus in p^ which will 
tend to be peaked in the regions in between conducting 
layers. We note that periodicity and inversion symmetry 
imply that p{z) — p{—z) = p{z + c). 

We give a layer index to the random potential Vj de- 
pending on whether the electron z co-ordinate lies in 
layer 1 or layer 2. There will then be two important 
types of disorder correlations (Vi(r')y,(r)) where i and 
j are layer indicies. Firstly, correlations within a single 
layer, (Fi (r') Fi (r)) , which contribute to the in-plane elas- 
tic scattering rate, and appear in the disorder-averaged 
single layer Green's functions. Secondly, correlations be- 
tween layers (Vi(r')V2(r)), that are relevant for interlayer 
conductivity. The corresponding expressions are (using 
the symmetry properties of Pz): 



(V^i(r')^i(r)) 
{Vi{v')V2{v)) 



L2 



d'v" 



dzpz{z) 



xl^„,p(r' - r", z)M,„p(r - r", z), (12) 



- / Sr" 



dzpz{z) 



xV^mp(r' - r",c - z)Mmp(r - r", -z), 

(13) 



C. Elastic scattering times in a layered metal 





In order to be clearer as to what wc mean by an inter- 
layer Cooperon, we now discuss the nature of the disor- 
der potential and the associated elastic scattering times. 
We shall assume that disorder correlations can be repre- 
sented as delta correlated, i.e., 



(y,(r')F,(r)) = a,^.<5(r-r'), 
where i, j = 1, 2. It is useful to rewrite 



(14) 



t/2 = 



^3 



1 



m*Tij 2Tij nNs 



where iV, 



27rh' 



is the two-dimensional density of states, 



not including spin degeneracy, and m* is the effective 
mass. Thus there are two distinct elastic scattering times 
associated with the disorder, th = T22 is the scattering 
time for in-plane elastic scattering, whilst T12 = T21 is 
the scattering time corresponding to scattering correla- 
tions between adjacent layers. If the impurities are all 
equidistant between the layers, i.e., p{z) has a maximum 
at z — ±c/2, then we will have rn ~ T12 , whereas if 
the impurities are located in the conducting layers, we 
expect Til < ^12. 

We note that previous theoretical treatments of AMRO 
have not considered the possible role of disorder correla- 
tions between neighbouring layers which could also mod- 
ify the semi-classical interlayer conductivity.— 



A, 2 



A, 2 



FIG. 1: Diagrams that give rise to the "interlayer Cooperon", 
described by Eq. (|17p . Neighbouring layers are indicated as 
1 and 2. A and R denote advanced and retarded Greens 
functions, respectively. 



D. Interlayer Cooperon in the absence of magnetic 

field 



The conductivity diagrams corresponding to the situ- 
ation we consider are as shown in Fig. [21 Following the 



FIG. 2: Diagrams for the quantum correction to the interlayer 
conductivity. Note that C — U12C is related to the interlayer 
Cooperon. This is a diagrammatic representation of Eq. (I16p . 



treatment of the two-dimensional problem by Rammer,— 
or Datta)^ we can write the correction to the interlayer 
conductivity from the Cooperon (in zero magnetic field) 



as 



(15) 6a^^ 



e-^cti 



de 



d^n / d^r2 / d'rg / d^u 



df 

x{Gf(ri,r3)Gf(r4,r2)C7f2Ci2(r3,r4) 
xG^(r2,r3)G^(r4,ri) + Gf(r4,ri)Gf(r2,r3) 
xC/i'2C2i(r3,r4)G2«(ri,r3)G2«(r4,r2)} , (16) 

where G = UI2C is the Cooperon [Fig.[T]. The diagram- 
matic expansion of the Cooperon in Fig. [1] implies the 
following equation for the interlayer Cooperon: 

Ci2(r,r') = <5(r-r') + y'd2r"jf2(r,r")Gi2(r",r'),(17) 

where 

jf2(r,r') = i7i22Gf(r,r')G^(r,r'), (18) 

is the Cooperon insertion. We can Taylor expand 
Gi2(r",r) and obtain^- the following generalization of 
the standard equation for the Cooperon in the experi- 
mentally relevant limit that the in-plane elastic scatter- 
ing time. Til, is much less than the inelastic scattering 
time, Til <C T0: 



I- - DaV^X Gi2(r,r') = —5{v - r') 

iTa J Til 



(19) 



which has the same form as the equation for the intralayer 
Cooperon, except the diffusion constant is 



Til 



•Jp'll 



T12 2Ti2 



(20) 



and the scattering rate which destroys the interlayer 
Cooperon is 

11 1 1 

TaiT) Til T12 T^iT) 



(21) 



In contrast, the first two terms do not appear for the 
intralayer Cooperon. For a general impurity potential 
Til < ■''12 £md so the effective dephasing rate will have 
a non-zero value at zero temperature. The structure of 
this expression has some similarities to what occurs with 
spin-orbit scattering^^ or intervalley scatteringi^ 



E. Interlayer Cooperon in the presence of a 
magnetic field 

In the presence of a magnetic field parallel to the layers 
we note that the Green's functions in layers 1 and 2 are 
related by a gauge factor— 



Gi(r,r')=exp(^— B||.(r-r')JG2(r,r'), 

and that the Green's functions can be written as a gauge 
dependent phase multiplying a gauge invariant piece: 

G(r,r') =exp f^J\^.dl\ Go(r,r'), 

where A^ is the vector potential associated with the per- 
pendicular magnetic field B^ and then 



Ji^,(r,r') ^ Ji^,'"(r,r')exp(^B||.(r-r') 



2h 



2ie, 



xexp(— (r-r')-A^(r)j. (22) 

This allows us to generalize Eq. p^ . as follows. We 
choose the Landau gauge to determine the perpendic- 
ular field, make a gauge transformation to eliminate 
the Bii dependent term, and end up with an equation 
for G;2(r,r') = exp(^B|| • (r - r')) Gi2(r, r') which is 
(when we approximate J12' by its zero field value) 



Da y 



2ie 



C[,ir',r) = —S{r~r'). 

Til 



(23) 



tC,0 



Note that the approximation of J12' by its zero field value 
implies that for high enough fields, our expressions will 
not be applicable. Now, Gi2(r3,r4) ~ Gi2(r3,r3), since 
the Cooperon is dominated by loops of time-reversed 
paths, and in a perpendicular magnetic field B±, the 
solution of equations of the form of Eq. ([23]) is well 
knowni^Si^ For our parameters this is: 



^I2(r3,r3) = ^f 



1 



2^^ ^0 4^'^l^l^^'r(" 



1) + 111 



(24) 



with rtmax — ^s/^^r^ where £ = vpTu is the in-plane 
mean free path. 



To calculate the conductivity, we need to take into ac- 
count the phase acquired by the propagators due to the 
vector potential differing in each layer. Thus, considering 
the first term in the conductivity (since the second term 
is its complex conjugate) and switching to momentum 
space propagators the expression simplifies to: 
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TThL"^ 
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G«(ki)G^(k2) 
B,i 



/ rfr3 j dr4Gi2(r3, r4)e'(''i+''^)-('-^--^) 



(25) 



Making use of Eq. ([24|) we obtain the expression: 



5a1 = 



2ceH\Uf2 feB 



Trh 
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Re[F(B||)] 
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(26) 



In Appendix |^ we show that for ecBn ^ hkp (the ex- 
perimentally relevant limit) 



where 



Re[F(B||)] 



/(^) 



1 



f TBnce 



2TO*Tfi 

, 1 + X 

' (1 + ^)3/2- 

^ for a; ^ 1 and f{x) 



(27) 

(28) 
i for 



Note that f{x) 
a:> 1. 

We can rewrite Eq. (|26p in a similar form to the mag- 
netoresistance in a thin film^^ using the properties of the 
digamma function Tp(x),— which is related to the sum we 
have obtained via 



tp{x + n + 1) — ip{x) — 2^ 



1 



771—0 
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Hence we obtain 
5a1{B^) = -ai 
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(29) 



with Bo = -^, and 



2ceHi 



(30) 



where Bp is the field scale defined in Eq. ^. 

The magnetic field dependence is of the same form as 
that for the weak localization correction to intralayer con- 
ductivity, Eq. (O, with B^{T) replaced by the magnetic 
field scale 



Ba{T) 



AeDaTairy 



(31) 



where Tq is given by Eq. (|2ip . Consequently, we note 
that unhke the phase breaking field associated with the 
intralayer conductivity that this quantity will be non-zero 
even at zero temperature. In the absence of a perpendic- 
ular magnetic field the interlayer Cooperon correction to 
conductivity reduces to 



SafiB^ = 0,r) = criln 



Tl2 

raiT) 



(32) 



Thus, the interlayer resistivity increases logarithmically 
as the temperature decreases. 

We now compare these quantum corrections to the 
semi-classical conductivity which describes conventional 
AMRO. The latter is given byi^ 



J„ (7 tan 61)2 



(Tconv(iJ) =Cro > -j— — 7 

^^ 1 + (nCJoTii 



n— — 00 



COS( 



(33) 



where 7 = ckp, Jn{z) is the Bessel function of order n, 
ujQ = -S^-, and the conductivity at zero field and at 6* = 
is 



o-Q 






(34) 



We can estimate the relative size of the quantum correc- 
tion to the conventional term as 



ctq 



Trkpl 



1 



(35) 



for fields such that eB\\ct/h <C 1. This indicates that 
the effects we discuss here are most likely to be seen in 
"dirty" samples (i.e., short mean- free path) and with a 
small Fermi surface. 



Peak in the vicinity of ^ = 7r/2 



Now, rimax depends on polar angle - as9 



B, 







and hence Ur, 



CO. Using the asymptotic expansion 



of "0(2;) at large x, we obtain (in the limit that ^ 
^ » 1) 
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-CTi <^ In 



Bo 
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Bl 



'^^Bl\ 



(36) 



This correction will have two different dependences on 
parallel field arising from / (|i?||c£). At fields much less 
than Bp — h/ecl it will be independent of parallel field, 
whilst for fields larger than Bp it takes the form 



Sa 



c 



CTO Bp 



irkpf. B, 



Bl 



^"'t)-24BM 



(37) 



In the limit that 9 



we can also compare with the 



limit from Uconv which can be obtained by standard sad- 
dle point methods applied to an integral representation 



ofEq. 



In the limit that 7 tan 9^1 
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O-Q 



Bp 

B^i 



1 + 2 cxp 



Wo Til cos 9 



sin(27tan( 



If we compare the two contributions in the limit that 
6* -^ -I, for fields larger than Bp, then we find that the 
size of the quantum correction relative to the semiclassi- 
cal magnetoconductivity is the same as for the zero field 
correction and has the same temperature dependence 



Sa^ 



1 



■In 



nkpi \Ta{T) 



T12 



(38) 



In order to compare with experiment, it is useful to 
separate the conductivity as cr(B) — CTback(B) -I- Act(B), 
where aback (B) = o'conv(B) + Sa^'{B± — 0), since it is 
only the negative magnetoresistance for increasing Bj_ 
that will be visible. Defining pback = 1/cback, we have 



Pzz 
Pback 



1 



Act 



c 



fback 



which implies a decrease in the resistivity, which grows as 
B±^ increases, i.e. as 9 deviates from ^. In consequence, 
there will be a peak above the background resistance at 



IV. DISCUSSION 

The correction to the resistivity at = ^ we found in 
the previous section is much smaller in magnitude than 
the conventional term, but since conventional AMRO are 
featureless in this region, this quantum contribution can 
still be visible above the background. The Cooperon 
should lead to a small peak in the resistivity in the vicin- 
ity of 61 = ^. The values of 9 for which this should be 
true are those such that both asymptotic formulae hold, 
i.e., cos 61 ^ 1/(c/cf) and Bcos9 = B± ^ Ba- 



A. Numerical Results 

We now evaluate the resistivity correction obtained in 
the previous section, Eq. (I29p . for parameters appropri- 
ate for a-(ET)2KHg(SCN)4 (Refs. llli). Fig. |1 shows 
the calculated temperature dependence and Fig. [5] the 
angular dependence. We calculate the dependence of the 
interlayer resistivity on temperature, field, and field di- 
rection. We choose kp = 1.4 nm~^, ep = 40 meV, th = 2 
ps, T12 = I.OIth and c ~ 20 A. We choose t^ — 66 ps, 
which gives Tq = 50 ps, a reasonable value based on the 
discussion following Eq. ([7]). From Eq. ([7]) we see that 
T0 ~ 66 ps at r ~ 1.2 K. These choices correspond to 
fields of Bo '- 0.2 T and Ba '- 0.002 T. The width of the 
peak is approximately A9 = 0.5° and the height of the 
peak is about 0.5% of the background resistance at all 
fields. 




FIG. 3: Interlayer resisitivity as a function of temperature 
for several different magnetic fields perpendicular to the lay- 
ers. In the absence of the field the resistivity increases log- 
arithmically with temperature due to weak localization ef- 
fects. This feature disappears as the magnetic field becomes 
comparable to the field Ba- The resistivity is normalized to 
its value at 0.2 K. The parameters used in these plots are 
kp = 1.4 nm^^, c = 2nm, rn = 2 ps, ri2 = l.Olrn, ep = 40 
meV, and t^{T)/tii — 50{K/T)^ . For these parameter values 
Ba{T = 0) =0.002 T. 



0.998 



Sr 0.996 



0.994 



0.992 



0.99 




FIG. 4: Interlayer resisitivity as a function of the angle 6, 
between the magnetic field and the normal to the layers. The 
peak at 90 degrees is due to weak localisation effects. Note 
that it is only about 0.5 percent of the total resistivity. The 
resistivity is normalized to the sum of the semi-classical resis- 
tivity and the zero field quantum correction. As the magnetic 
field increases from 2 to 20 T, the size of the peak decreases 
and its width increases. The parameters used in these plots 



are kp = 1.4 nm 



— 2nm, Til = 2 ps, ri2 — l.Olrn 



ep = 40 meV, and Ta — 50 ps. 



In Fig. Owe replot the data from Fig. |4]as a function 
of B± ~ Bcos{6). It is clear that as B increases, the 
peak becomes a function of B± only. We determined the 
width of the peak in two different ways: i) by choosing a 
fixed value of p{9)/p{9 — ^), and then determining the 
appropriate A9 as a function of B, and ii) by subtracting 
the (non-flat) background and finding the value of A9 
corresponding to the half-maximum of the peak. In case 
i) A9 r^ B~^ for essentially all values of B, whereas for 
case ii) A9 - B-^ for B > Sp ~ 6 T. 

The results in Figs. and [HI are qualitatively in accord 
with recent experiments on the resistivity in /3"-(BEDT- 
TTF)4[H30Cr(C204)3]CH2Cl2, in which a peak that de- 
pends only on B cos{9) for small 6 and a peak width 
A9 ~ 1/B were observed. ^^ Additionally, features in un- 
published data^ on a-[(ET)i_:r(BETS):r]2KHg(SCN)4 
with X ~ 0.03, appear to be consistent with the picture 
presented here. 
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FIG. 5: Resistivity data from Fig. |4] plotted as a function 
of the magnitude of the magnetic field parallel to the layers 
Bcos9. 



B. Implications for experiment 

Inelastic scattering plays an important role in deter- 
mining when a resistivity peak will be observed. A rea- 
sonable temperature dependence of the inelastic scatter- 
ing rate implies that at high temperatures, the inelastic 
scattering time is much less than the elastic scattering 
time r^ <C Til. At a temperature where the two become 
comparable, one might expect to see the first hints of a 
resistance peak which narrows and grows with decreas- 
ing temperature until the peak width and height satu- 



rate. We expect there will be a saturation when I/tq 

reaches its low temperature limit ^ — — which is 

generally non-zero. The visibility of the peak is largest 
when Ta 3> Til, which will be true when T12 is close in 
magnitude to th. 

In concert with the magnetic field dependent signa- 
tures discussed above, the temperature and magnetic 
field dependence of such a peak in the resistivity should 
allow it to be distinguished from coherence peaks, whose 
width and height should have much weaker dependence 
on temperature and magnetic fieldf^^ As a consequence, 
these localization peaks also offer an opportunity to de- 
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FIG. 6: Width of the peak in the resistivity about 6 — 90° 
as a function of field for the data in Fig. U) The width is 
determined either by a fixed drop in the resistance, or by sub- 
tracting the background. Fits to both procedures are shown. 



magnetic field and the inelastic scattering time, which 
should allow it to be distinguished from coherence peaks 
seen in very low disorder layered metals. This potentially 
allows for the extraction of information about inelastic 
scattering in layered metals, and gives another tool with 
which to learn about the physics of these systems. 
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termine the inelastic scattering rate, independent of reg- 
ular AMRO which allow a determination of the elastic 
scattering rate, and quantum oscillations which allow a 
determination of the total scattering rate through the 
Dingle temperature. 

Some questions that we have not addressed here but 
provide interesting avenues for further exploration regard 
the effects of anisotropy and the coherence of interlayer 
transport. Our calculation was for an isotropic Fermi 
surface with isotropic elastic scattering - in general one 
can expect Fermi surface properties to vary on different 
parts of the Fermi surface, which would likely imply that 
the nature of the peak for parallel fields will also de- 
pend on the orientation of the magnetic field within the 
plane. This might allow the possibility of investigating 
the inelastic scattering on different parts of the Fermi 
surface. Secondly, our derivation assumed weakly inco- 
herent interlayer transport. However, it appears likely 
to us that similar results should also hold for coherent 
interlayer transport, and so the presence or absence of 
this weak localization peak is not necessarily evidence 
for or against coherent interlayer transport. The visibil- 
ity of the peak will be stronger in samples with relatively 
smaller values of kpl, corresponding to smaller values of 
T, but the criterion usually used to argue coherence of in- 
terlayer transport is comparing the hopping amplitude, 
t±, with the scattering rate -. 



V. CONCLUSIONS 

In conclusion, we have calculated the effect of the in- 
terlayer Cooperon on AMRO in quasi-two dimensional 
layered metals. The interlayer Cooperon can give rise to 
a peak in the magnetoresistance for fields parallel to the 
layers that behaves similarly to conventional weak local- 
ization. The features of the peak depend sensitively on 



APPENDIX A: CALCULATION OF F(B ) 

To calculate ^ (By ) we need to perform an integral over 
four in-plane disorder averaged Green's functions 



xG^(ki 



2h 



B|| G-4 -ki 



(2^)2 [E-Ek 
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2h 
1 



B 



IT] J E - ek-p-irj 

(Al) 



with 77 = h/2Tii, ek = l;^, and ek-p = efe - 

J^l3cos(l>+^, with P = \ecB\\. Rewriting Eq. (|XT1) 
as an integral over energy and performing the energy in- 
tegral using contour integration, we get 
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The second term in the denominators is much smaller 
than the third term if /? ^ 2hkp. This corresponds to 
B\\ < 4:hkp/{ec) ~ 0(100 T), which is always satisfied 



for physically realistic parameters. We then have 



where 



1 /-^^ 
2^ Jo ^ + ^. 



cos^ 



(A4) 
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Evaluating this by contour integration gives Eqn. 
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